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1. 


Formulas  and  tables  are  provided  making  it  possible  for  users  of  the 
truncated  sequential  probability  ratio  tests  of  MIL-STD  781C  to  make 
acceptance  or  rejection  decisions  when  the  test  procedure  is  prematurely 
truncated. 


2 .  Introduction 

Users  of  the  truncated  sequential  probability  ratio  tests  given  in 
MIL-STD  781C  [2]  are  typically  concerned  with  sampling  a  lot  of  items  each 
of  which  is  assumed  to  have  identical  exponential  life  distributions.  A 
number  n  of  items  are  kept  on  test;  as  soon  as  one  fails,  it  is  replaced 
by  another.  Let  0  denote  the  mean  time  between  failures.  If  6  ■  9q 
the  lot  is  considered  acceptable;  if  9  »  0^ ,  9j  <  8q ,  the  lot  is  consid¬ 
ered  unacceptable.  Let  N(t)  denote  the  number  of  failures  by  time  t 
where,  throughout,  t  is  measured  in  units  of  9^.  The  truncated  MIL-STD 
781C  tests  provide  two  boundaries,  an  upper  boundary  A(t)  and  a  lower 
boundary  B(t)  and  a  truncation  time  T.  If  by  time  T,  N(t)  crosses 
A(t)  before  it  crosses  B(t),  then  the  hypothesis  9*9^  is  accepted;  if 
N(t)  crosses  B(t)  before  it  crosses  A(t),  or  neither  boundary  is 
crossed  by  time  t  ■  T,  then  the  hypothesis  9  ■  0q  is  accepted.  A 
producer's  risk  a  (probability  of  rejecting  9q  when  it  is  true)  and  a 
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consumer’s  risk  0  (probability  of  accepting  6{1  when  Is  true)  are 

provided  for  each  of  the  eight  tests  in  MIL-STD  78 1C. 

Consider  the  following  problem.*’  A  lot  is  being  tested  using  one  of 
the  eight  tests  of  MIL-STD  781C.  At  some  time  t  (t  <  T)  it  becomes 
necessary  to  make  a  determination  about  the  lot;  l.e. ,  to  either  accept  or 
reject  the  lot.  If  N(s)  has  crossed  either  A(s)  or  B(s)  at  some 
s  £  t,  then  the  decision  is  dictated  by  the  MIL-STD  78 1C  test  being  used. 
However,  if  such  is  not  the  case,  the  user  of  MIL-STD  781C  has  no  basis  for 
making  a  decision.  We  present,  here,  a  procedure  that  permits  prematurely 
terminated  MIL-STD  781C  sequential  testing  and  construct  tables  and 
decision  rules  that  can  be  used  to  supplement  MIL-STD  781C  for  uae  in  such 
situations.  Our  calculations  are  based  on  the  method  presented  by  Kao, 

Kao,  and  Mogg  [1]. 


3.  Prematurely  Terminated  Sequential  Teats  (PTST) 

For  t  <  T,  let  p^  denote  the  probability,  when  0  is  the 

true  mean  time  between  failures,  that  N(t)  ■  k  and  N(s)  has  not  crossed 
either  A(s)  or  B(s)  for  s  £  t.  Each  truncated  MIL-STD  781C  test  has 
associated  sets  of  values  {u^}  and  {v^};  u^  is  the  solution  of 
A(u^)  ■  i  and  v^  is  the  solution  to  BCv^  -  i  where  the  i’s  are 
Integers.  These  numbers,  called  reject  and  accept  numbers  in  MIL-STD  781C, 
are  listed  with  each  MIL-STD  781C  test.  They  depend  on  the  particular  test 
used.  Using  the  method  in  [1],  it  is  possible  to  calculate,  for  each  test, 


*The  problem  was  presented  by  William  Wallace  of  the  Naval  Electronics 
System  Command  at  the  Quality  Control  Workshop  held  in  Washington,  D.C.  on 
on  March  7-8,  1983. 
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associated  sets  of  values  (A^(d/9^)}  and  (B^(9/9j)};  the  range  of 
indices,  1,  of  {Aj(0/9^)}  are  the  same  as  for  {u^}  and  the  range  of 
indices,  i,  of  {B^(0/0^)>  are  the  same  as  for  {v^}.  The  interpretation 
of  B1(9/0p  is  that  it  is  the  probability,  when  0  is  the  true  mean  time 
between  failures,  of  N(t)  crossing  B(t)  at  t  *  not  having  previ¬ 
ously  crossed  either  boundary  (N(t)  can  cross  B(t)  only  at  points 
(v^).  A^O/9^),  also,  has  a  boundary  crossing  probability  interpretation; 

however,  for  our  purposes  this  does  not  have  to  be  made  explicit.  The 
equation  for  ^9/9^)  is  given  by 


where  p  equals  one  plus  the  greatest  Integer  less  than  A(0)  and  q'(t) 
is  the  greatest  integer  less  than  B(t).  Let  ct  be  a  critical  number  at 
time  t;  i.e.,  if  it  should  be  necessary  at  time  t  to  accept  or  reject 
9  ■  9n,  then  accept  9*9,.  if  N(t)  <  c  .  Then, 


(2)  P{Accepting  0q  by  time  t ] ©>  -  P{Accepting  9^  prior  to  time  t J  0> 

+  P{Accepting  0q  at  time  t  j  Q} 


q’(t)  ct 

-  I  B  <8/0  )  +  l 

j-0  J  k-[q'(t)]+l 


Pk,t<e/0l) 


t 


where.  In  the  case  t  <  Vq,  the  lower  limit  of  summation,  [q ’ ( t ) ]  +  1,  is 
interpreted  as  zero. 

Given  c£,  a  test  based  on  the  termination  time  t  and  the  original 
truncated  MIL-STD  781C  test  is  defined.  We  refer  to  this  as  a  prematurely 
terminated  sequential  test  (PTST).  The  producer's  risk  aCc^)  is  given  by 

(3)  a(ct)  -  1  -  P{Accepting  0Q  by  time  t|0Q} 


and  the  consumer's  risk  0(ct)  is  given  by 

(4)  0(ct)  -  P{Accepting  0Q  by  time  tjOj}  . 


We  stress  that  for  (3)  and  (4)  to  be  valid,  the  termination  time  t 
must  be  independent  of  N(s),  s  <  t.  For  otherwise,  a  value  of  t  could 
be  selected  to  favor  one  or  the  other  hypothesis. 

If,  for  each  t,  cfc  is  chosen  such  that  a(ct)  £  a,  the  producer's 
risk  of  the  truncated  MIL-STD  781C  in  use,  then  the  producer's  risk  for  the 
overall  PTST  procedure,  allowing  for  any  premature  termination  time  t, 
which  is  Independent  of  N(s),  s  £  t,  to  be  less  than  or  equal  to  a. 

Exact  equality,  if  desired,  can  be  achieved  by  a  suitable  use  of  randomiza¬ 
tion.  In  the  same  way  the  consumer's  risk  can  be  controlled.  However, 


both  risks  cannot  be  simultaneously  guaranteed. 


There  are  two  other  problems  that  may  be  of  interest.  The  user  may  be 


concerned  with  the  conditional  probability  that  he  accepts  0g  given  that 
he  has  been  forced  to  make  a  decision  at  time  t  and  N(s),  s  <  t,  has  not 
crossed  a  boundary.  This  conditional  probability  is  simply 


I 

k-(q'(t)]+l 


Pk.t^/V7 


I 

B(t)<k<A(t) 


Pk,t<e/61> 


The  user  may  also  be  Interested  in  more  flexibility  in  choosing, 
initially,  a  truncation  time.  The  PTST  provides  the  means  for  constructing 
additional  tests  with  a  greater  variety  of  truncation  times. 


4.  Tables  of  Risks  for  PTST 

For  each  of  the  eight  truncated  sequential  tests  of  MIL- STD  781C  we 
calculate  {A^CB/Q^)} ,  {A^(l)},  and  {B^(l)>.  For  values  of 

t  ■  1/5  T,  ...»  4/5  T  we  tabulate  p^  t(0j/0j)  (i  *  0,1)  for  the  rele¬ 
vant  values  of  k.  We  define  c(t,a)  to  be  the  smallest  value  of  c£ 
such  that  a(cfc)  <  a.  We  find  c(t,a)  and  calculate  a(t)  -  a( c(t,a)) 
and  p(t)  ■  0(c(t,a)). 

If  t  t  1/5  T,  ... ,  T,  then  one  can  approximate  by  using  the  largest 
multiple  of  T/5  less  than  t  for  his  termination  time  or  the  formulas 
(l)-(4)  can  be  used  to  derive  the  PTST  for  that  particular  time  t. 

5.  Example 

Suppose  a  shipment  of  electronic  equipment  is  to  be  tested.  The  life 
of  each  piece  of  equipment  is  assumed  to  be  exponential  with  mean  p.  It 


would  be  acceptable  If  p  -  -  100  hours  and  unacceptable  if  p  -  Pj  •  67 

hours.  It  Is  decided  to  let  n  ■  10  so  that  0g  *  10  and  8j  ■  6.7. 

Test  Plan  IC,  which  has  an  a  ■  0  *  .10  and  a  T  *  49.50  measured  In  units 
of  6.7  hours,  is  to  be  used.  Suppose  the  test  is  prematurely  terminated  at 
t  ■  20  (134  hours).  Since  this  is  not  one  of  our  values  of  T  in  our 
tabulation,  we  round  down  to  2/5(49.5))  *  19.8.  From  the  tables  the 
critical  number  at  19.8  is  18.  Thus,  the  test  used  is:  Accept  0g  if 
prior  to  T  *  19.8,  N(s)  crosses  B(s)  before  crossing  A(s)  or  if 
N(19.8))  <  18.  The  value  of  ct(19.8)  -  .085.  The  value  of  p(19.8)  -  .388 
If  it  is  relevant  to  consider  t  as  a  premature  termination  necessi¬ 
tating  an  early  decision  given  that  no  decision  has  been  reached;  the 
conditional  probability  that  the  hypothesis  0  ■  9^  is  rejected  (when 
9  -  90)  ia 


21 


21 


Jl9  -  -0726 


k-il 


When  0  -  9,  the  conditional  probability  that  the  hypothesis  8  *  0q  is 
accepted  is 


18 


21 


j  K.19.S ■  -72”  • 

c-11  k»ll  ’ 


These  conditional  probabilities  can  be  interpreted  as  the  risks  associated 
with  only  premature  decisions  (excluding  those  cases  where  decisions  may 
have  been  reached  using  the  MIL-STD  procedure). 


If,  in  fact,  0  -  10,  the  conditional  probability  that  the  hypothesis 


0  •  0q  ■  10  is  rejected  given  N(s)  has  not  crossed  either  A(s)  or  B(s) 
prior  to  t  ■  19.8  is 


21 


I 

k-19 


21 


Pk  19  3<1.5)/  I 
k’19*8  k-11 


pk,19.8 


(1.5) 


.0726 


and  if,  in  fact  0  -  6.7,  the  conditional  probability  that  the  hypothesis 
0  -  0q  «  10  is  accepted  given  N(s)  has  not  crossed  either  A(s)  or 
B(s)  prior  to  t  «  19.8  is 


18 


I 

k-11 


21 


pk,19.8(1)/ J 


k-11 


pk,19.8 


(1) 


.7299  . 
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10 

0.0248 

0.0525 

11 

0.0052 

0.0221 

9 

0.0359 

0.0048 

10 

0.0518 

0.0140 

11 

0.0448 

0.0242 

12 

0.0287 

0.0311 

13 

0.0133 

0.0288 

14 

0.0027 

0.0116 

t,.l) 

a(t) 

P(t) 

4 

.066 

.590 

7 

.093 

.369 

11 

.103* 

.271 

14 

.114* 

.193 

Although  a  »  .10  is  the  nominal  producer's  risk  In  MIL-STD  781C  IIIC 
our  method  of  calculation  of  the  MIL  STD  producer's  risk  yields  an  a 
which  exceeds  our  tabulated  values. 


MIL-STD  78 1C  TEST  IVC 


a  ■ 

.20 

0  - 

.20 

0  - 

eo 

91  * 

.5  e0 

1 

A1(2) 

BjU) 

A^l) 

Bt(l) 

• 

0 

0.2466 

0.0608 

1 

0.1732 

0.0428 

2 

0.0432 

0.1158 

0.1217 

0.0284 

" — 

3 

0.0513 

0.0781 

0.1451 

0.0192 

- 

4 

0.0380 

0.0523 

0.1077 

0.0129 

•m, 

5 

0.0276 

0.0778 

6 

0.0194 

0.0548 

y  .y: 

7 

0.0137 

0.0388 

t 

k 

Pk,t<2> 

>k,t(1> 

1.95 

0 

0.3776 

0.1426 

1 

0.3678 

0.2777 

2 

0.1560 

0.2356 

3.90 

1 

0.1996 

0.0569 

2 

0.2403 

0.1371 

3 

0.1371 

0.1563 

4 

0.0245 

0.0559 

5.84 

3 

0.1364 

0.0587 

4 

0.1062 

0.0915 

5 

0.0365 

0.0629 

7.79 

4 

0.0688 

0.0224 

5 

0.0728 

0.0474 

6 

0.0362 

0.0471 

7 

0.0028 

0.0074 

t 

c(t,.2) 

a(t) 

8(t) 

1.95 

2 

.099 

.656 

3.90 

3 

.176 

.411 

5.84 

5 

.185 

.345 

7.79 

7 

.206* 

.276 

Although  a  *  .20  13  the  nominal  producer's  risk  In  MIL-STD  78 1C  IVC, 
our  method  of  calculation  of  the  MIL  STD  producer's  risk  yields  an  a 
which  exceeds  out  tabulated  value. 
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MIL- STD  78 1C  TEST 

VC 

a 

m 

.10 

8  -  .10 

0 

m 

®0 

0X  -  .33  8q 

Ai(3) 

Bi(3) 

AjU)  B1(l) 

0.2865 

0.0235 

0.2066 

0.0170 

0.0149 

0.1384 

0.0919  0.0113 

0.0275 

0.0911 

0.1689  0.0074 

0.0200 

0.1230 

0.0138 

0.0845 

0.0094 

0.0575 

k 

Pk,t( 

3) 

Pk,t(1> 

0 

0.5016 

0.1262 

1 

0.3461 

0.2612 

2 

0.1103 

0.2498 

1 

0.3145 

0.0597 

2 

0.2329 

0.1327 

3 

0.0902 

0.1541 

4 

0.0072 

0.0370 

2 

0.1831 

0.0262 

3 

0.1576 

0.0678 

4 

0.0708 

0.0903 

5 

0.0125 

0.0483 

3 

0.1048 

0.0113 

4 

0.1042 

0.0338 

5 

0.0535 

0.0521 

6 

0.0136 

0.0396 

MIL-STD  781C  TEST  VIC 


Ul  W  IO  H 


MIL-STD  78 1C  TEST  VI IC 


A1(1.5) 


0.0398 

0.0425 

0.0367 


B1(1.5) 

0.1225 

0.1140 

0.0948 

0.0763 


A1(l) 


0.0893 

0.0957 

0.0826 


BjU) 

0.0429 

0.0399 

0.0332 

0.0267 


'k,t(l*5) 

0.4039 

0.2567 

0.3662 

0.3491 

0.1660 

0.2374 

0.0140 

0.0299 

0.1631 

0.0659 

0.2958 

0.1792 

0.2682 

0.2437 

0.1475 

0.2010 

0.0238 

0.0488 

0.1383 

0.0533 

0.2310 

0.1334 

0.2124 

0.1840 

0.1244 

0.1617 

0.0278 

0.0542 

0.1041 

0.0382 

0.1792 

0.0986 

0.1699 

0.1403 

0.1046 

0.1296 

c(  t , .3) 


MIL-STD  78 1C  TEST  VIIIC 
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